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Abstract. We investigate the cosmological behavior in a universe governed by time asym¬ 
metric extensions of general relativity, which is a novel modified gravity based on the addition 
of new, tinre-asymmetric, terms on the Hamiltonian framework, in a way that the algebra of 
constraints and local physics remain unchanged. Nevertheless, at cosmological scales these 
new terms can have significant effects that can alter the universe evolution, both at early 
and late times, and the freedom in the choice of the involved modification function makes 
the scenario able to produce a huge class of cosmological behaviors. For basic ansatzes of 
modification, we perform a detailed dynamical analysis, extracting the stable late-time so¬ 
lutions. Amongst others, we find that the universe can result in dark-energy dominated, 
accelerating solutions, even in the absence of an explicit cosmological constant, in which the 
dark energy can be quintessence-like, phantonr-like, or behave as an effective cosmological 
constant. Moreover, it can result to matter-domination, or to a Big Rip, or experience the se¬ 
quence from matter to dark energy domination. Additionally, in the case of closed curvature, 
the universe may experience a cosmological bounce or turnaround, or even cyclic behavior. 
Finally, these scenarios can easily satisfy the observational and phenomenological require¬ 
ments. Hence, time asymmetric cosmology can be a good candidate for the description of 
the universe. 
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1 Introduction 

The standard model of cosmology includes two accelerated phases of expansion, at early and 
late times respectively. Such a behavior cannot be obtained within the standard paradigm 
of physics, namely in the framework of general relativity and Standard Model of particles. 
Hence, additional degrees of freedom should be included in the picture. If these extra degrees 
of freedom are attributed to new, exotic ingredients of the universe content, then concerning 
late times one has the concept of dark energy (for reviews see [1, 2]) and concerning early 
times the concept of inflaton field(s) (for reviews see [3, 4]). On the other hand, if the extra 
degrees of freedom are of gravitational origin, then one obtains the paradigm of modified 
gravity (see [5, 6] and references therein). The latter approach has the additional motivation 
of improving the UltraViolet behavior of gravity and alleviating the difficulties towards its 
quantization [7, 8]. Note that there are not strict boundaries between the above approaches, 
since one can partially or completely transform from one to the other, or construct theories 
where both extensions are imposed. 

In the usual approach to gravitational modification one adds higher-order corrections to 
the Einstein-Hilbert action, like in F(R ) gravity [9-13], in Gauss-Bonnet and f(G ) gravity [14, 
15], in Lovelock gravity [16, 17], in Weyl gravity [18, 19], in Horava-Lifshitz gravity [20-22], 
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in Galileon modifications [23-26], in nonlinear massive gravity [27-30] etc. A different class 
of gravitational modifications arise when one starts from the equivalent torsional formulation 
of gravity and add higher-order correction, like in f(T) gravity [31-34], in f(T,Tc ) gravity 
[35-37], etc. 

Recently, a new class of modified gravity was proposed [38]. In particular, working in 
the Hamiltonian framework the authors constructed a theory that breaks the time reversal 
invariance of general relativity. Although the algebra of constraints and local physics are 
unchanged, new terms appear at cosmological scales, that can alter the universe evolution, 
both at early and late times. 

In the present work we are interesting in investigating in detail the cosmological impli¬ 
cations of the above time asymmetric extensions of general relativity. In order to achieve 
this independently of the initial conditions and the specific universe evolution, we apply the 
dynamical systems method [39, 40] which allows us to extract the global behavior of the 
scenario, bypassing the complexity of the involved equations. Indeed, due to the freedom in 
choosing the relevant extra modification function, the capabilities of the scenario are found to 
be huge. The plan of the work is the following: In section 2 we present the time asymmetric 
extension of general relativity and we apply it in a cosmological framework. In section 3 
we perform a detailed dynamical analysis, extracting the stable late time solutions and the 
corresponding observables, and in section 4 we discuss their physical implications. Lastly, 
section 5 is devoted to the conclusions. 

2 Time asymmetric extensions of general relativity and cosmology 

Let us briefly review the time asymmetric extension of general relativity [38]. In a first 
subsection we present the gravitational model itself, while in a second subsection we apply 
it in a cosmological framework. 

2.1 Time asymmetric extension of general relativity 

In this formulation one starts with the Hamiltonian form of general relativity with a cosmo¬ 
logical constant [41] 

S GR = JdtJ {tt ab g ab - NU ADM - N a V a } , (2.1) 

where 

u ADM = ~y/g{R - 2A) + (vr a \ab - ^ 2 ) + M* = 0 (2.2) 

is the usual Hamiltonian constraint. In the above expressions g a b is the spatial metric, with 
7 r ab its canonical momenta and 7r = g a b ft ab the corresponding trace, while N and N a are the 
usual lapse and shift functions. In this formalism, the Hamiltonian constraint (2.2), along 
with the diffeomorphism constraint 

V a = D b Tr b a + V* = 0, (2.3) 

form a first class algebra, where the terms and T>^ correspond to the matter content and 

D a is the covariant derivative. Obviously, the above expressions respect the time reversal 
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symmetry 


t —^ —t 


(2.4a) 

(2.4b) 

(2.4c) 


9ab ^ 9ab 


7r 


ab 


—7r 


ab 


In order to acquire well defined cosmological evolution equations one must use a gauge 
fixing, and it proves convenient to use the “constant mean curvature gauge condition” (CMC) 

[38] 

vr - y/g < 7T >= 0, (2.5) 

where < • • • > denotes the spatial average of a density p defined through < p >= (J s p) / (J 2 yd)), 
with V = / v yjg the spatial volume. The CMC condition (2.5) is a gauge fixing of the Hamil¬ 
tonian constraint (2.2), and thus they form a second class system. However, note that the 
CMC condition (2.5) and the diffeomorphism constraint (2.3) form a system of four first 
class constraints [42-44], as it is the case for the Hamiltonian constraint along with the dif¬ 
feomorphism constraint. One can show that, restricting to constraints that are local in g a b 
and 7 r ab , there are no other pairs of systems of four first class constraints that one is the 
gauge fixing of the other, however one has the freedom to add a term linear in 7r to the 
Hamiltonian constraint [44], This new term 7r/L, with L the length-scale where this term 
becomes significant, breaks the time reversal symmetry (2.4a)-(2.4c), and this feature gave 
to the obtained gravitational modification the name “time asymmetric extension of general 
relativity”. One can extend the above extra, time-asymmetric, term of the Hamiltonian 
constraint, by assuming that the length-scale in which it becomes important is driven by a 
function of spatially averaged quantities, such as the spatial volume V. Hence, in summary, 
one can extend (2.2) to a modified Hamiltonian constraint of the form [38] 

n new = -^ y /g(R- 2A) + ^ (Aa;, - ^tt 2 ) + / {V)ir + 77* = 0, (2.6) 

where f(V ) is an arbitrary function of V. 

The above modification of the Hamiltonian constraint gives rise to a novel class of grav¬ 
itational modifications. The new term leaves the constraint algebra and the local physical 
degrees of freedom unchanged [38]. The only complexity comes from the fact that it affects 
the propagation of chiral fermions, since the left-handed spacetime connection D a ^a does 
depend on 7r ab . In order to handle this issue, one introduces the Ashtekar geometry [45], 
alongside the usual spacetime geometry characterized by the spacetime metric g^ u . Thus, 
although the gravitational effects and the propagation of photons are governed by the conven¬ 
tional spacetime geometry, the propagation of chiral fermions is determined by the Ashtekar 
geometry which contains all the information of time irreversible behavior. Nevertheless, since 
in this work we are interested in the late-time background cosmological evolution, in which 
the matter sector is effectively described by a perfect fluid, and where radiation (a part of 
which is composed by chiral fermions) is negligible, in the following we do not discuss the 
above issue in more details. Hence, the time asymmetric modified gravity that we focus in 
this work is characterized by the action 

S-Jdtf {t r ab g ab - Nn new - N a V a } , (2.7) 

where is given by (2.6). 
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2.2 Cosmological application of time asymmetric gravity 


Let us now apply the time asymmetric extension of general relativity in a cosmological 
framework. In particular, we focus on a Friedmann-Robertson-Walker (FRW) spacetime 
metric of the form 


ds 2 


—dt 2 + a(t ) 2 


dr 2 

1 — kr 2 



( 2 . 8 ) 


where a(t) is the scale factor, k = —1,0,1 for spatially open, flat or close geometry respec¬ 
tively, and with dfl 2 the two-dimensional sphere line element. Note that tinre-asymmetric 
extension of general relativity singles out a specific 3 + 1 decomposition, selected by the 
constant mean curvature gauge condition, and moreover it introduces a dependence on the 
spatial slices volume, and thus the spacetime must be spatially compact. This is indeed the 
case in the above cosmological metric, where k refers to positive, negative or zero constant 
spatial curvature. In particular, all of these cases are consistent with a spatially compact 
topology, with k = +1 corresponding to spheres, k = 0 to tori, while for k = —1 the infinite 
number of compact manifolds with constant negative curvature are classified by Thurston 
[46]. Inserting the above metric in the total action S + S m , with S given by (2.7) and S m 
the matter action, and performing the variation in the ADM formalism, we easily obtain the 
Friedmann equations as [38] 


k 8 ttG 


H z + ^ = 
a z 


pm + f(V(a)) z 


(2.9) 


H — ~ 2 = -47T G(p m + Pm) + af(V(a)) d ^V( a ^ ; (2.10) 

a z oa 

where H = a/a is the Hubble parameter, V(a) oc a 3 is the spatial volume, and G is the 
gravitational constant. Additionally, we have considered the matter action S m to correspond 
to a perfect fluid with energy density p m and pressure p m respectively. We stress here that 
in action (2.7) we do not include an explicit cosmological constant, since our goal is exactly 
to investigate whether the universe acceleration can arise solely from a general modification 
term f(V) (which definitely in the specific case f(V) = const, gives rise to an effective 
cosmological constant). 

Defining for convenience g(a) = ^f(V(a)), the above modified Friedmann equations 
become 


k 87 tG 


H z + — = 


'Pm + 


G 2 g(a ) 


H-\ - - 4 *G( Pm+ p m ) + 

a z a a z 

and thus the modification is included in the arbitrary function g(a). Furthermore, 

rewrite the Friedmann equations (2.11),(2.12) in the usual form 

H 2 = y (■ P™ + Phz) 


H = — — (Pm +p m + PDE + PDE ) > 


( 2 . 11 ) 

( 2 . 12 ) 

we can 


(2.13a) 

(2.13b) 


if we define the energy density and pressure of the effective dark energy sector as 

3 G g(a) 2 


Pde — 


87t a 2 
G 


PDE = 


9 (+ 2 9( a )9'( a ) 

o ' 


(2.14) 

(2.15) 
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i.e. attributing the dark energy sector to the new terms that time asymmetric gravity brings 
to the Friedmann equations. In this case, the dark energy equation-of-state parameter be¬ 
comes: 


wde 


Pde 

Pde 


1 

3 


2 ag'jay 

9(a) . ’ 


(2.16) 


In summary, the modified gravity at hand is determined by the arbitrary function 
/(17(a)). Hence, according to the choice of /(17(a)) one obtains distinct classes of cosmolog¬ 
ical models. 


3 Late-time cosmology 

In this section we are interested in investigating in detail the late-time cosmology of the time 
asymmetric extension of general relativity. Since the gravitational modification is determined 
by the function f(V), we will choose two basic ansatzes, namely the power law and the 
exponential one. In particular, we will consider 

• Model I: f(V) = giV m , which implies that the auxiliary function g(a) becomes g(a) = 
% (^) > w ith p = 3 m + 1, with g\ a constant and p a parameter, and where ao is a 
constant which can be set to 1 for convenience. 

• Model II: f(V ) = g 2 e xv , which implies that g(a) = < 72 §e Aa3 / a o, with g 2 a constant, A a 
parameter, and with ao a constant which can be set to 1. 

In order to study the cosmological behavior in a general way, independently of the 
initial conditions and the specific universe evolution, we will apply the dynamical systems 
method, which allows to extract the global features of a cosmological scenario [47-54]. In 
this procedure, one first transforms the involved cosmological equations into an autonomous 
system and then he extract its critical points. Hence, perturbing linearly around these 
critical points, and expressing the perturbations in terms of a perturbation matrix, allows to 
determine the type and stability of each critical point by examining the eigenvalues of this 
matrix. 


3.1 Model I: f(V) = g l V m 

In the case where f(V ) = g\V m , i.e. when g{a) = |kaP (with p = 3m +1), with g\ a constant 
and p a parameter, the Friedmann equations (2.11),(2.12) become 

2 & _ 8nG 2 2p—2 /q i \ 

H -| 2 — q Prn + 5l ° > (3.1) 

CL O 

H —\ = -4 irG(p m + p m ) + (p - l)gfa 2p ~ 2 , (3.2) 

a z 

and thus the effective dark energy (2.14) and pressure (2.15) respectively become 


Pde = 


3 gf 

8ttG C 


2p—2 


PDE = ~^Q a2P 2 ( l + 2 P)’ 


(3.3) 

(3.4) 
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and hence (2.16) leads to 


1 

vjde = -^(1 + 2 P>- 


(3.5) 


Additionally, we can define the “total” equation-of-state parameter 


_ 2 H _ 8irGa 2 w m p m - g\ (2p + l)a 2p + k 

tot 3 H 2 8ir Ga 2 p m + 3 gfa 2 P — 3k 

and the deceleration parameter as 


H 1 + 3 w to t 
q ~ 1 H 2 ~ 2 ’ 


(3.6) 


(3.7) 


with w m = pm/Pm the matter equation of state. In the following we assume the usual energy 
conditions, which lead to 0 < w m < 1. Finally, note that for g\ = 0 we re-obtain standard 
general relativity. 


3.1.1 Zero or negative curvature 

In the case k = 0, —1 as auxiliary variables it proves convenient to use the various density 
parameters, namely 


- 1 /,: — 


k n _ 87TG Prn n _ gla 2p 2 

a 2 H 2, m 3 H 2 ' llDE ~ H 2 


and thus the first Friedmann equation (3.1) gives rise to the constraint 


+ Qde = 1 - 


(3.8) 


(3.9) 


Using the above auxiliary variables we can write the cosmological equations in the au¬ 
tonomous form 

— 7 "“ = — [2 P^IdE + ( 3u) m + 1 )(Qde + - 1)], (3.10a) 

arj 

, DE = —II DE [2 p{^DE — 1) + (3 Wm + 1 ){^DE + Hfc — 1)] , (3.10b) 

dr] 

where we have used the constraint (3.9) in order to eliminate and thus reduce the 
system to dimension two. In these equations, as usual, we define the logarithmic time 
r] = lno. Hence, the above autonomous system is defined on the compact phase space 
{(fife, £Ide) ■ Hfc > 0 ,^de > 0, flfc + Qde < 1} 1 - Finally, using the auxiliary variables (3.8) 
we can express the deceleration parameter (3.7) as 

q = 2 [! + 3 Wm^m - (2 P + 1 )£Ide - Hfe] • (3.11) 

The scenario of Model I, namely f(V) = giV m , with zero or negative curvature, ad¬ 
mits three physical critical points, corresponding to expanding universe (H > 0), which are 
displayed in Table 1 along with their existence conditions. In the same Table we include 

x The interest of defining compact phase spaces is that then the flow has well-defined past and future 
attractors, and this facilitates the drawing of global results for the cosmological scenario [47-54]. 


- 6 - 











C.P. 


Qde 

q 

Existence 

Eigenvalues 

Stability 

Pi 

0 

0 

' 6 w m + l 

2 

always 

3ui m + l,2p + 3w m + 1 

saddle for p < 3w ™ +1 

unstable for 0 < w m < 1 ,p > — 3m ™ +1 

P 2 

1 

0 

0 

always 

2 p, —3w m - 1 

stable for p < 0 
saddle for p > 0 

Pi 

0 

1 

-p 

always 

-2 p, -2 p - 3 w m - 1 

unstable for p < — 3w ™ +1 
saddle for — 3w ™ +1 < p < 0 
stable for p > 0 


Table 1. The physical critical points of the system (3.10) of time asymmetric cosmology of Model I: 
f(V) = g\V m , with zero or negative curvature, and their existence and stability conditions. We have 
assumed 0 < w m < 1. 


the eigenvalues of the involved perturbation matrix, and thus the corresponding stability 
conditions. Finally, for completeness, we also include the values of the deceleration param¬ 
eter, calculated through (3.11). Note that the solution associated to P 3 for p / 1 is the 

1 

power-law form a(t ) = [(1 — p ) (ci + ait)] , while for p = 1 it is just the de Sitter solution 
a(t) = c±e 9lt , with ci and ai integration constants. 

In summary, the scenario at hand admits two stable late-time critical points, namely 
Pi for p < 0 and P 3 for p > 0. 


3.1.2 Positive curvature 


In the case k = +1, that is for positive curvature, it is not guaranteed that the Hubble 
parameter does not change sign during the evolution. This implies that the H-normalization 
that we used in the previous open and flat case is not a good choice for creating compact 
variables, since when H crosses zero the dynamical variables would diverge, and moreover 
when H change sign our “time” variable 77 = In a would change flow. Thus, in the present 
k = +1 case, it is consistent to introduce the auxiliary variables (similarly to the variables 
introduced in section VI of [39], and in sections 3.3 and 5.3 of [55]) as: 


©fc 


1 

^PV 2 ' 


Q 0 


H 

V 


, &r 


87T Gp m 
3V 2 ’ 


©DE 


g 2 a 2 ^ 

V 2 


(3.12) 


where V = ■s/H 2 + a 2 , and which are finite even if H crosses zero. Therefore, the first 
Friedmann equation (3.1) leads to the constraint 


©m + &DE = 1- 

Additionally, from the definition of D it follows 

©fc + Qo = 1 ) 

while the curvature parameter is expressed as 

1 


(3.13) 

(3.14) 


Hz. — 


1 -Ql 


a 2 H 2 Q 2 


Using the above auxiliary variables we can re-write the cosmological equations as 

(3.15a) 


= - (l — Qq) [2 P@DE + (3 Wm + !)(©£>£ — 1)] , 


d@DE 

dr 


= — Qo(2p + 3 w m + 1)(@de ~ 1 )&DE, 


(3.15b) 
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where we have used the constraints (3.13) and (3.14) in order to eliminate 0 m and and 
therefore reduce the system to dimension two. In the above dynamical system, we have 
introduced the consistent “time” variable r through dr = Ddt, which indeed satisfies the 
necessary requirement that it is monotonic even if H change sign. The above autonomous 
system is defined on the compact phase space {( Qo,®de ) : — 1 < Qo < 1,0 < @de < 1 }. 
Finally, using the auxiliary variables (3.12) we can express the deceleration parameter (3.7) 
as 

q = xvvj [1 + 3w m (l — &de ) — (2 p + 1)©de] ■ (3.16) 

2C^q 

As we can observe the system (3.15) is symmetric under the transformation 

(r, Qo, ®de) —> (—t, — Qo, ®de)- (3-17) 

Thus, it is sufficient to discuss the behavior in one part of the phase space, that is in t > 
0 , Qo > 0 ,&de > 0 , and then obtain the dynamics on the other part from (3.17). For 
example, if a point with coordinates (Qq, @* de ), Qq > 0,®* DE > 0 is a future attractor as 
r —>• +oo, then its partner point (—Qo, ®*de) v ^ a (3-17) is a past attractor as r —>• —oo, and 
vice versa. Furthermore, we mention that the function 

M = V 9 ^ ~J~ = ~(^ w m + l)Q 0 Af, (3.18) 

X C^q cLt 

is a monotonic function in the regions Qo < 0 and Qo > 0 for ®de 7^ 1- The points having 
Qo > 0 correspond to expansion, while those having Qo < 0 correspond to contraction. As 
we will see in the following, the system (3.15) admits a fixed point with Qo = 0 if we assume 
p > 0,0 < w m < 1. However, since for p < 0,0 < w m < 1 there are not equilibrium points 
with Qo = 0, it follows that M acts as a monotonic function in the interior of the phase 
space. As a consequence, for p < 0 there can be no periodic orbits in the interior of the phase 
space and global results can be drawn [39]. Additionally, from the definition of M it follows 
that either Qq —>• 1 or 0 de —>• 1 asymptotically. 

Note that the system (3.15) allows for an easy analytical elaboration, leading to 

Til-(“l gCi 

Qo( a ) = 1 - a 2p+3wm+l + e ci 5 ®DE(a) — 1 - a 2 P+3Wrn+ i + gC1 ) (3-!9) 

with ci and ai integration constants. 

The scenario at hand admits five physical critical points which are displayed in Table 
2 along with their existence conditions. In the same Table we include the eigenvalues of 
the corresponding perturbation matrix, and the resulting stability conditions. Finally, we 
also include the values of the deceleration parameter, calculated through (3.16). Note that 
the points P4 and P5 have the time reversal behavior of Pq and Pj respectively, due to the 
symmetry (3.17). Additionally, there exist orbits connecting P 4 and P$ with Pq and Pj, 
which implies that Qo can indeed become zero, i.e. H = 0 , during the evolution (recall that 
the points having Qo > 0 are expanding while those having Qo < 0 are contracting). Lastly, 
the system admits a static solution, namely P%, which always behaves as a saddle point. 

In summary, the scenario of Model I, namely f(V ) = g\V m , with positive curvature, 
admits three stable late-tinre critical points, namely the expanding solution P 5 for p > 0 , 
the contracting solution Pq for p > —(3 w m + l )/2 and the contracting solution P7 for p < 
~(3 w m + l)/2. 
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C.P. 

Qo 

P>DE 

q 

Existence 

Eigenvalues 

Stability 

Pi 

l 

0 

2 

always 

3 w m + 1, 2p + 3 w m + 1 

saddle for p < 3u, ” +1 

unstable for p > — 3m ” +1 

Po 

l 

1 

-p 

always 

—2 p, -2 p - 3 w m ~ 1 

unstable for p < — 3w ™ +1 
saddle for — 3w ™ +1 < p < 0 







stable for p > 0 

p e 

-l 

0 

iWm+i 

2 

always 

— (3Wm + 1). -(2 P + 3 Wm + 1) 

saddle for p < 3u, ” +1 

stable for p > 3w ™ +1 

P-1 

-l 

1 

-p 

always 

2 P, 2 P + 3Wm + 1 

stable for p < 3w ™ +1 

saddle for — 3w ™ +1 < p < 0 







unstable for p > 0 

Po 

0 

3Wm+l 

undefined 

p > 0 

-V , p( 1 + 3w m), vM 1 + 3w m ) 

saddle 

-*8 

2p+3tu m +l 


Table 2. The physical critical points of the system (3.15) of time asymmetric cosmology of Model I: 
f(V) = g\V m 1 with positive curvature, and their existence and stability conditions. 


We close this paragraph with some comments on the auxiliary variables choice. The 
advantage of using the variable Qo versus using the variable fl*., used in paragraph 3.1.1, is 
that for closed models the variable would not keep track of the H -sign changes, due to 
the quadratic dependence on Qq , however these changes may have important cosmological 
consequences. As we saw, choosing Q o instead of fl*. allows us to to obtain novel features, 
such as expanding solutions, contracting partners, transition from contracting to expanding 
cosmologies and vice versa, as well as static solutions (for instance a static solution, where 
H = 0, would obviously not be seen using 17-normalization). These differences, arising from 
the possible H -sign change in closed models, forbids a unified description of all cases. Such 
a difference between closed and open/flat geometries, and the implied necessary different 
normalization, was first observed in [39, 56, 57] despite the fact that closed FRW had been 
previously studied in [48, 58, 59]. 


3.2 Model II: f(V) = g 2 e xv 

In the case where f(V) = g 2 e xv , i.e. when g{a) = g 2 ^e Xa3 , with g 2 a constant and A a 
parameter, the Friedmann equations (2.11),(2.12) become 


2 k 8irGp m 2 „ 2 Aa 3 

ti + — = — t, -f- g 2 e 


a* 


k 


H - - = -4:TrG(p m + p m ) + 3a s g2\e 


3 „2 \ n 2\a 3 


and thus 


(3.20) 

(3.21) 


vj£)E = —1 — 2Aa 3 . 

Additionally, the “total” equation-of-state parameter reads 

_ 2 H 8TrGa 2 w m p m — 3g\a l e 1Xa3 (1 + 2Aa 3 ) + k 

tot 3 H 2 8TrGa 2 p m + 3 g 2 a 2 e 2Xa3 — 3k 

while the deceleration parameter writes as 

1 -{- 3U!/o/, 

Finally, note that for g 2 = 0 we re-obtain standard general relativity. 


(3.22) 


(3.23) 


(3.24) 
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3.2.1 Zero or negative curvature 

In the case k = 0, —1, we introduce the density parameters compact auxiliary variables 




k n _ 8 tt Gp m 0 _ g%e 2a3x 

a 2 H 2 ’ m 3iL 2 ’ H 2 


and thus the first Friedmann equation (3.20) gives rise to the constraint 


^k + + HflE = I- 


(3.25) 


(3.26) 


In order to be able to close the system we need one more auxiliary parameter. Since the 
corresponding choice proves to be different according to the sign of A, we will examine the 
two cases separately. 


• A > 0 

In this case we define the additional auxiliary variable 

A a 3 

T = - 

1 + Aa 3 


(3.27) 


Since by construction 0 < T < 1 (since A > 0), we can define 


^L = H(l-T)-\ (3.28) 

which implies fj = ^Aa(f) 3 + ln[a(t)j (modulo an additive constant), and thus fj —>• —oo 
as a —> 0 and fj —>• 00 as a —>• 00 . Hence, using the auxiliary variables (3.25) and (3.27) 
we can re-write the cosmological equations in their autonomous form, namely 


dT _ 
dfj 
dklk 
dfj 

d&DE 

dfj 


3T(1 — T) 2 , (3.29a) 

= 3(T — 1 )(w m + 1 )^k(^DE + Alfc — 1) — 2fife [T(3£Ide + — 1) — Q k + 1], 

(3.29b) 

= 3 (T — 1 ){w m + 1 )£Ide(Qde + Alfc — 1) + 2 Qde [^k — T(3 Qde + A4 — 3)], 

(3.29c) 


where we have used the constraint (3.26) in order to eliminate Clearly, the above 
system is defined on the { (T. Qde) ■ 0 < T < 1, > 0, kloE > 0, + &de < 1} 

part of the phase space, where we have included the two boundaries T = 0 and T = 1. 
Furthermore, note that the invariant subset boundary T = 1 corresponds to the asymp¬ 
totic future, while the invariant subset boundary T = 0 is associated asymptotically to 
the (classical) initial state. Therefore, in this formalism all the fixed points are located 
at T = 0 and T = 1 [60]. Lastly, using the auxiliary variables (3.25), (3.27) we can 
express the deceleration parameter (3.24) as 



1 T 3 


1 + T 
1 - T 


&DE 



(3.30) 


10 











Label 


Qde 

T 

<? 

Existence 

Eigenvalues 

Stability 

Qi 

0 

0 

0 

3 w m + 1 

2 

always 

3,3(1 + w m ), 1 + 3 w m 

unstable 

Q2 

1 

0 

0 

0 

always 

3, —3(1 + Wm), 2 

saddle 

Q3 

0 

1 

0 

-1 

always 

3, —3(1 4- w m ), —2 

saddle 

Qi 

0 

0 

1 

arbitrary 

always 

6 , 0,0 

nonhyperbolic, 

behaves as unstable 

Qh 

^kc 

0 

1 

arbitrary 

Qkc G ( 0 , 1 ] 

6 , 0,0 

nonhyperbolic, 

behaves as unstable 

Q% 

0 

1 

1 

—00 

always 

- 6 ,- 6,0 

nonhyperbolic, 

behaves as stable 


Table 3. The physical critical points of the system (3.29) of time asymmetric cosmology of Model II: 
/( V) = < 72 e AV , with zero or negative curvature and A > 0, and their existence and stability conditions. 
We assume 0 < w m < 1. 


Thus, for the critical points having T = 0, the expression (3.30) is well-defined and 
gives q |t=o = [1 + 3w m Q m — 3Q.de ~ 14] /2, while for the critical points having T = 
1 ,Qde > 0 we obtain q —>• —oo as T —>• l~ since Q m ,Q.DE and 14 are bounded. On 
the other hand, for the critical points having T = 1,Qde = 0, q is arbitrary. 

The scenario at hand admits five physical critical points, and one curve of critical 
points (namely Q 5 ), which are summarized in Table 3 along with their existence condi¬ 
tions. In the same Table we include the eigenvalues of the corresponding perturbation 
matrix, and the resulting stability conditions. Finally, we also include the values of 
the deceleration parameter, calculated through (3.30). We mention that for the three 
nonhyperbolic critical points, the linear analysis is not adequate to determine their 
stability, and therefore the stability conditions have been extracted applying the center 
manifold method [61]. The corresponding investigation is performed in Appendix A.l. 


• A < 0 

In this case we define the additional auxiliary variable 

U - -2^. (3.31) 

Since 0 < T\ < 1 (since A < 0), we can define 

^ = H(1-T 1 )~ 1 , (3.32) 

which implies 77 = —^Ao(t) 3 +ln[a(t)] (modulo an additive constant), and thus 77 —>• —00 
as a —>• 0 and fj —>• 00 as a —>• 00 . Hence, using the auxiliary variables (3.25) and (3.31) 
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Label 

^k 

&DE 

Ti 

q 

Existence 

Eigenvalues 

Stability 

Qt 

0 

0 

0 

2 

always 

3,3(1 +W m ),l +3 Wm 

unstable 

Q& 

1 

0 

0 

0 

always 

3, —3(1 + Wm), 2 

saddle 

Q9 

0 

1 

0 

-1 

always 

3, -3(1 + Wm), -2 

saddle 

Q 10 

0 

0 

1 

arbitrary 

always 

- 6 , 0,0 

nonhyperbolic, 
behaves as stable 

Q11 

^kc 

0 

1 

arbitrary 

ftfcc € ( 0 , 1 ] 

— 6 , 0,0 nonhyperbolic, 

behaves as saddle for Ofc c 7 ^ 1 . 
stable for = 1 . 

Ql 2 

0 

1 

1 

+00 

always 

6 , 6,0 

nonhyperbolic, 
behaves as unstable 


Table 4. The physical critical points of the system (3.33) of time asymmetric cosmology of Model II: 
f(V ) = < 72 e AV , with zero or negative curvature and A < 0, and their existence and stability conditions. 
We assume 0 < w m < 1. 


we can re-write the cosmological equations in autonomous form as 


-p- = 3Ti(l -Ti) 2 , 
dr/ 

^ = 3Q k [(Tr - l)w m 


(3.33a) 

1 ] iS^DE + — 1 ) + [Ti(9£Ide ~ 1 ) — 2 ] + (T\ + 2 ) 0 |, 

(3.33b) 


d&DE 

df) 


^de {3 (flpE — 1) [Ti(wm + 3) - w m - 1] + (Ti - l)(3u; m + l)fl fc } , (3.33c) 


where we have used the constraint (3.26) in order to eliminate Clearly, the above 
system is defined on the {(Ti, Q k ,^DE) ■ 0 < Ti < 1, Q k > 0, flpE > 0, 0, k + £Ide < 1} 
part of the phase space, where we have attached the two boundaries T\ = 0 and T\ = 
1. Finally, note that in terms of the auxiliary variables (3.25),(3.31) the deceleration 
parameter (3.24) is given by 


q 


1 

2 


1T 3w m Q m 3 


O) 


^DE ~ ^k 


(3.34) 


Thus, for the critical points having T\ = 0, the expression (3.34) is well-defined and 
gives q \t 1= o = [1 + 3 vu m £l m — 3f Ide ~ flfc] /2, while for the critical points having Ti = 
1 ,Qde > 0 it follows that q —>• +oo as Ti —>• 1” since and Q k are bounded. 

On the other hand, for the critical points having Ti = 1, Qde = 0, q is arbitrary. 

The scenario at hand admits five physical critical points, and one curve of critical 
points (namely Q n), which are summarized in Table 4 along with their existence condi¬ 
tions. In the same Table we present the eigenvalues of the involved perturbation matrix 
and the corresponding stability conditions. Finally, we also include the values of the 
deceleration parameter, calculated through (3.34). Concerning the three nonhyperbolic 
critical points the stability conditions have been extracted applying the center manifold 
method [61]. The corresponding investigation is performed in Appendix A.2. 

In summary, the scenario of Model II, namely f(V) = g 2 & XV , with zero or negative 
curvature admits the following stable late-time solutions: For A > 0 the expanding solution 
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Qq (nonhyperbolic, but with stable center manifold). For A < 0 the line of expanding 
solutions Qio (nonhyperbolic, but with stable center manifold) and a point on the line Qu 
with flfcc = 1 (nonhyperbolic, but with stable center manifold) which represent a curvature 
dominated solution. 


3.2.2 Positive curvature 

In this case we introduce the density parameters compact auxiliary variables (similarly to 
the variables introduced in section VI of [39], and in sections 3.3 and 5.3 of [55]): 


©fc 


1 


Qo 


H 

V' 


Qr 


8 irGp m 
W 2 ’ 


&DE 


g ¥- 2 “ 3A 

v 2 ' 


(3.35) 


with V 


y/H 2 + a 2 , and thus the first Friedmann equation (3.20) gives rise to the constraint 


©m + &de = 1, (3.36) 

and as before we have the restriction 

©fc + Qo = l- (3.37) 


In order to be able to close the system we need one extra auxiliary parameter. Since the cor¬ 
responding choice is different for different signs of A, we will examine the two cases separately. 


• A > 0 

In this case we define the additional auxiliary variable 


1 + Aa 3 

Since by construction D > 0 and 0 < T < 1 (since A > 0), we can define 


dr 

dt 


V{1 -T)- 1 , 


(3.38) 


(3.39) 


which implies r = J * ^ + - d( (modulo an additive constant). Since T>{ 1 — 

T) _1 > 0, f is a monotonic function of t. Hence, using the auxiliary variables (3.35) 
and (3.38) we can express the cosmological equations as 


— = 3Q 0 T(1-T) 2 , (3.40a) 

^ = \ (1 - Qo) [-3(T - 1 ){w m + 1 )(G de - 1) + T(6 @de - 2) + 2] , (3.40b) 

~ DE = 3Qo(@DE - 1 )&DE [ T(w m - 1) - w m - 1] , (3.40c) 

dr 

where we have used the constraint (3.36) in order to eliminate © m and (3.37) to elim¬ 
inate ©a,. The above system is defined on the Qq,Qde) ■ 0 < T < 1, — 1 < Qo < 

l,0<©DE<lf part of the phase space, where we have included the two boundaries 
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T = 0 and T = 1. For Qo > 0, i.e. for expanding cosmologies, the invariant sub¬ 
set boundary T = 1 corresponds to the asymptotic future, while the invariant subset 
boundary T = 0 is associated asymptotically to the (classical) initial state. However, 
for contracting models (Qo < 0) the roles of the invariant sets T = 1 and T = 0 are 
reversed in time. This arises from the fact that the function 

N = N' = 3(1 — T)QqN , (3.41) 

is a monotonic function in the region Qo < 0,0 < T < 1, where N is monotonically 
increasing, and in the region Qo >0,0 <T <1, where N is monotonically decreasing, 
and thus it follows that either T —>• 1 or T —>• 0 asymptotically. Hence, in this formalism 
all the fixed points are located at T = 0 and T = 1 [60]. 

Furthermore, the system (3.40) is invariant under the symmetry 

(f ,T, Q 0 , ®de ) —>■ (—t, T, —Qo, ®de) ■ (3.42) 


Thus, it is sufficient to discuss the behavior in one part of the phase space, for instance in 
T >0, T > 0, Qo > 0, ®de > 0, while the dynamics on the other part is being obtained 
from (3.42). For example, if a point with coordinates (T*, Qjj, @* de ), Qq > 0,®* DE > 0 
is a future attractor as f ^ Too, then, its partner point (T*,—Qq,Q* de ) via (3.42) is 
a past attractor as f —>• — oo, and vice versa. 

Finally, using the auxiliary variables (3.35) and (3.38) we can express the deceleration 
parameter (3.24) as 


Q = 


2 Ql 


1 + 2>Wm{l - @DE ) — 3 


1 + T 
1 -T 


®DE 


(3.43) 


Thus, for the critical points having T = 0, the expression (3.43) is well-defined and 
gives q |t=o = [1 + 3rc m (l — &de) ~ 3 @de\, whereas for the critical points having 

T = 1, @de 7^ 0, it is implied that q —>• —oo as T —>• 1~. On the other hand, for the 
critical points having T = 1,Qde = 0, q is arbitrary. 

Furthermore, we note that the function 


M= 1 9 ^f , M' = —{2>w m + 1)(1 — T)QqM, (3.44) 

1 v 0 

is a monotonic function in the regions Q o < 0, 0 < T < 1 and Qo>0,0<T<l for 
®de / 1- The points having Qo > 0 are expanding, and those having Qo < 0 are 
contracting. Thus, from the definition of M it follows that either Q(j— >-lor0£>£;—>-1 
asymptotically. 

The scenario at hand admits four physical critical points, and one curve of critical 
points (namely Qi6, which is the straight line joining the points (0,0,1) and (1,0,1), 
with the left point not included) corresponding to expanding cosmologies, which are 
summarized in Table 5 along with their existence conditions. In the same Table we 
include the eigenvalues of the involved perturbation matrix, and the corresponding 
stability conditions. Moreover, we also include the values of the deceleration parameter, 
calculated through (3.43). Each of the above critical points have contracting partners 
via the discrete symmetry (3.42), which are displayed in Table 6. 
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Label 

Qo 

Qde 

T 

<? 

Existence 

Eigenvalues 

Stability 

Ql3 

1 

0 

0 

3Wm + i 

2 

always 

3,3(1 + w m ), 1 + 3 w m 

unstable 

Q 14 

1 

1 

0 

-1 

always 

3, -3(1 + w m ), -2 

saddle 

Q 15 

1 

0 

1 

arbitrary 

always 

6,0,0 

nonhyperbolic, 
behaves as saddle 

Q 16 

Qoc 

0 

1 

arbitrary 

Qoc £ (o, i] 

6Q 0c ,0,0 

nonhyperbolic, 
behaves saddle 

Ql7 

1 

1 

1 

—oo 

always 

-6,-6,0 

nonhyperbolic, 
behaves as stable 


Table 5. The physical critical points of the system (3.40) of time asymmetric cosmology of Model II 
corresponding to expanding cosmologies: /( V) = g 2 & XV , with positive curvature and A > 0, and their 
existence and stability conditions. We assume 0 < w m < 1. 


Label 

Qo 

Qde 

T 

q 

Existence 

Eigenvalues 

Stability 

Rl 3 

-1 

0 

0 

3ui m + I 

2 

always 

-3, -3(1 + W m ), -(1 + 3 Wm, 

) stable 

Ru 

-1 

1 

0 

-1 

always 

-3,3(1 + w m ), 2 

saddle 

Rl 5 

-1 

0 

1 

arbitrary 

always 

-6,0,0 

nonhyperbolic, 
behaves as saddle 

Rie 

Qo c 

0 

1 

arbitrary 

Qoc e (0, l] 

—6Qoc 0,0 

nonhyperbolic, 
behaves as saddle 

R\7 

-i 

1 

1 

—OO 

always 

6,6,0 

nonhyperbolic, 
behaves as unstable 


Table 6. The physical critical points of the system (3.40) of time asymmetric cosmology of Model 
II corresponding to contracting cosmologies: /( V) = (? 2 e AV , with positive curvature and A > 0, and 
their existence and stability conditions. We assume 0 < w m < 1. 


Lastly, there exists a line of static solutions, i.e. neither expanding nor contracting, 
given by 


S\ : (: T,Q 0 , &de) 



(T c - l)(3w m + 1) \ 

3 T c {w m - 1) - 3 (w m + 1 )J ' 


T c 6 [0,1]. 


(3.45) 


Imposing the physical condition 0 < w m < 1, it follows that the above line always 
satisfies the existence condition 0 < Qde < 1. The eigenvalues of the linearization 
around S± are 


0 V (3 w m + 1) (1 - T c ) [( w m - 1)T C -w m -l] [(2 w m + l)r| - (w m + 6)T C ~w m -l] 

('Wm - 1 )T C -W m - 1 

V 7 (3 w m + 1) (1 - T c ) \{w m - l)T c \{2w m + l)r| - (w m + 6)r c -w m -l] 

(■Wm - 1 )T C -w m -l 

and thus whenever this line exists these eigenvalues are always real. Since two of them 
have different sign the whole line behaves as saddle. 

We mention that in order to determine the stability of the six nonhyperbolic critical 
points (expanding and contracting ones) we apply the center manifold method [61], and 
the corresponding analysis is performed in Appendix A.3. 
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• A < 0 

In this case we define the additional auxiliary variable 


^ _ A « 3 
Tl 1 - Aa 3 ■ 


(3.46) 


Since 0 < T\ < 1 we define 

df 

— = V(1-T 1 )~\ (3.47) 

so that f is monotonic increasing for T>(1 — Ti) -1 > 0. Therefore, using the auxiliary 
variables (3.35) and (3.46) we can re-write the cosmological equations in autonomous 
form as 


— 1 = 3Q 0 Ti(l-T 1 ) 2 , (3.48a) 

^ ~ (1 - Ql) [3(Ti - 1 )(w m + 1)(@ DE - 1) + 2(3T 1 @ de + Tx - 1)] , (3.48b) 

\ DE = 3Qq(®de — l)®DE{Ti{w m + 3) - w m — 1), (3.48c) 


where we have used the constraint (3.36) in order to eliminate fl m . The above system 
is defined on the {(Ti, Qq, @de) ■ 0 < Ti < 1, — 1 < Qq < 1,0 < @de < 1} part of the 
phase space, where we have attached the invariant boundaries T\ = 0 and T\ = 1. 

Similarly to the previous section, the function 


N = -^-, N' = 3(1 — Ti)QqN , 
t — b 


(3.49) 


is a monotonic function in the region Qq < 0,0 < T\ < 1, where N is monotonically 
increasing, and in the region Qq > 0,0 < Ti < 1, where N is monotonically decreasing, 
and it follows that either Ti —> 1 or T\ —>• 0 asymptotically. 

Furthermore, the system (3.48) is invariant under the symmetry 


(f,Ti, Qq, @de) —>• (—t ,Ti,—Qq,®de) ■ (3.50) 


Thus, it is sufficient to discuss the behavior in one part of the phase space, that is in 
f > 0, Ti > 0, Qo > 0 ,Q ee > 0, while the dynamics on the other part is being obtained 
from (3.50). For example, if a point with coordinates (Ti*, Qq, @ ee ), Qq > 0, @* DE > 0 
is a future attractor as f — >• +oo, then, its partner point (T\*, —Qq, Q* DE ) via (3.50) is 
a past attractor as f —>• — oo, and vice versa. Additionally, in the invariant set Ti = 0 
we obtain the first integral 


(1 _ 

(1 - ® de Y® 1 ^ 


(3.51) 


with c an integration constant. 

Finally, note that in terms of the auxiliary variables (3.35),(3.46) the deceleration pa¬ 
rameter (3.24) is given as 


Q = 


1 

2Q§ 


1 + 3w m (l - ®de) ~ 3 


O) 


®DE 


(3.52) 
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Label 

Qo 

@DE 

T\ 

Q 

Existence 

Eigenvalues 

Stability 

Qis 

l 

0 

0 

iSiOm+l 

2 

always 

3,3(1 + w m ), 1 + 3w m 

unstable 

Q 19 

l 

1 

0 

-1 

always 

3, —3(1 + W m ), —2 

saddle 

Q 20 

l 

0 

1 

arbitrary 

always 

- 6 , 0,0 

nonhyperbolic, 
behaves as saddle 

Q 21 

Qoc 

0 

1 

arbitrary 

Qoc £ ( 0 , 1 ) 

— 6 Q 0c , 0,0 

nonhyperbolic, 
behaves as saddle 

Q 22 

l 

1 

1 

+OO 

always 

6 , 6,0 

nonhyperbolic, 
behaves as saddle 


Table 7. The physical critical points of the system (3.48) of time asymmetric cosmology of Model 
II representing expanding cosmologies: f(V) = g 2 £ XV , with positive curvature and A < 0, and their 
existence and stability conditions. We assume 0 < w m < 1. 


Hence, for the critical points having T\ = 0, the expression (3.52) is well-defined and 
leads to q |ti=o = yyw [1 + 3w m (l — @de) — 30_d_e], however for the critical points 

having Ti = 1,@de 7 ^ 0, it follows that q —» +oo as T\ —>• 1“. On the other hand, for 
the critical points having T\ = 1, ®de = 0, q is arbitrary. 

The scenario at hand admits four physical critical points, and one curve of critical 
points (namely Q 21 ) representing accelerating solutions (Qo > 0), which are displayed 
in Table 7 along with their existence conditions. In the same Table we present the 
eigenvalues of the corresponding perturbation matrix, and the resulting stability con¬ 
ditions. Finally, we also include the values of the deceleration parameter, calculated 
through (3.52). Each point/curve in Table 7 has a partner through the symmetry 
(3.50), representing a contracting cosmology (Qq < 0), which are displayed in Table 8. 

Lastly, the system (3.48) admits a line representing static solutions, i.e. neither ex¬ 
panding nor contracting, given by 


S 2 ■ {T\,Qo, &de) 



(3 w m + 1) (T c - 1) \ 

3 ( w m T c + 3T c -w m -l)J ’ 


(3.53) 


with eigenvalues 


0 VT C - 1a/(3 W m + 1) [(Wm + 3)r e -W m - 1] [(4 W m + 15)T e 2 - 5 (Wm + 2 )T C + Wm + 1] 

(w m + 3 )T C - w rn - 1 

\JT C -lyf (3 w m + 1 ) [(Wm + 3 )T C - w m - 1] [{4w m + 15)r c 2 - 5 (w m + 2 )T C + Wm + 1] 

{wm T 3)T C Wm 1 


This line exists and is physical, i.e. possessing 0 < @de < 1, for i) 0 < T c < 

5Wm — \/3w m (^SWm "h8)H-40H - 10 -I .1 • -I -I i ri -1 -t jji 

- x — 8 ^ m+3Q -, when the eigenvalues are real and 02 behaves as saddle, or 

..X 5w m — \/3w m 8)40“hi0 rjn ^ 1 i , i i 

11 ) --— o , o n - < +c S 4 , when two eigenvalues are purely imaginary. 


We mention that in order to determine the stability of the six nonhyperbolic critical 
points we apply the center manifold method [61], and the corresponding analysis is 
performed in Appendix A.4. 


In summary, the scenario of Model II, namely f(V) = g 2 & XV , with zero or negative 
curvature admits the following stable late-tinre solutions: For A > 0 the contracting solution 
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Label 

Qo 

Qde 

Ti 

q 

Existence 

Eigenvalues 

Stability 

Rl 8 

-1 

0 

0 

Su'm + l 

2 

always 

-3, -3(1 + Wm), “(I + 3 Wm) 

stable 

Rl 9 

-1 

1 

0 

-1 

always 

—3,3(1 + w m ), 2 

saddle 

R20 

-1 

0 

1 

arbitrary 

always 

6,0,0 

nonhyperbolic, 
behaves as saddle 

R21 

Qoc 

0 

1 

arbitrary 

Qoc G ( 0 , 1 ) 

6Q00 0,0 

nonhyperbolic, 
behaves as saddle 

R22 

-1 

1 

1 

+00 

always 

-6,-6,0 

nonhyperbolic, 
behaves as saddle 


Table 8. The physical critical points of the system (3.48) of time asymmetric cosmology of Model 
II representing contracting cosmologies: /( V) = (? 2 e AV , with positive curvature and A < 0, and their 
existence and stability conditions. We assume 0 < w m < 1. 


i?i 3 and the expanding solution Qn (nonhyperbolic but with stable center manifold), while 
for A < 0 the contracting solution Pig. 

4 Physical Implications 

Having performed a complete dynamical analysis of cosmological scenarios governed by time 
asymmetric extensions of general relativity , we can now proceed to the discussion of the 
physical implications. In particular, we focus on the stable late-time solutions, since these 
solutions can attract the universe at late times, independently of the specific initial conditions 
and the specific intermediate evolution. 

4.1 Model I: f(V ) = g x V m 

In the case where f(V) = giV m , i.e. when g(a) = | kaP (with p = 3m + 1), with g\ a constant 
and p a parameter, and with open or zero curvature, the scenario at hand exhibits the 
three critical points presented in Table 1. Point Pi corresponds to a dark-matter dominated 
universe (fl m = 1 ), that is non-accelerating (q > 0), however it is never stable and thus it 
cannot attract the universe at late times. Point P 2 corresponds to a universe governed by 
the curvature term (f= 1), which is neither accelerating nor decelerating (this is typical 
for curvature dominated solutions [62]). For p < 0 it can be stable, and thus it can attract 
the universe at late times (this is actually expected since for p < 0 the effective dark-energy 
term decreases faster than the curvature term, and hence the latter dominates). However, 
its observational features are disfavored by observations. Point P 3 is stable for p > 0 and 
thus it can be the stable late-time state of the universe. It corresponds to a dark-energy 
dominated, accelerating universe, where the dark-energy equation-of-state parameter (3.5), 
namely wde = — (1 + 2 p)/3, can lie either in the quintessence regime (for 0 < p < 1), or 
in the phantom one (for 1 < p), or either behave as an effective cosmological constant (for 
p = 1) giving rise to a de Sitter universe. These features make it a good candidate for the 
description of the universe, especially if 0.9 < p < 1.1, in which case —1.07 < wde < —0.93 
in agreement with observations [63]. We mention that the above behavior is obtained without 
the addition of an explicit cosmological constant term in the action, i.e. it is a pure effect 
of the novel, time-asymmetric theory. Finally, note that even when the effective dark energy 
lies in the phantom regime, the universe does not end in a Big Rip [64-67], or any other type 
of singularity [68], at finite time. 
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Figure 1. The phase-space behavior of time asymmetric cosmology of Model I: f(V) = giV m , with 
negative curvature, p = 0.9 (i.e. m = —0.033 ), and w m = 0. The shadowed region marks the 
unphysical part of the phase space. In this specific example the universe is led to the the dark-energy 
dominated, accelerating solution P 3 . 


In order to present the above behavior in a more transparent way, we evolve numerically 
the cosmological equations and in Fig. 1 we depict the corresponding phase-space behavior. 
The unphysical part of the phase space (in which the density parameters exceed one) is 
marked by the shadowed region. As we can see, in this specific example the universe results 
in the dark-energy dominated, accelerating solution P 3 . 

In the case of positive curvature, the model possesses five critical points, displayed in 
Table 2 . Amongst them, the points P 5 and Pq can be stable, and thus they can attract the 
universe at late times. P 5 corresponds to an accelerating, dark-energy dominated universe 
(f lp>E = &de = 1 since for this point T> —> H in (3.12)), in which the dark-energy equation- 
of-state parameter can lie either in the quintessence or in the phantom regime, or behave 
like an effective cosmological constant. Hence, it can be a good candidate for the description 
of the universe. On the other hand, P$ corresponds to a matter dominated, contracting 
solution, and as we mentioned before it could be an attractor too, but it does not describe 
accurately the universe at late times. Moreover, point P 7 has the reverse dynamical behavior 
of P 5 due to (3.17), i.e it corresponds to a contracting (Qo < 0), dark-energy dominated 
universe (£Ide = @de = 1), and it can be an attractor too. Similarly, point P 4 presents 
the time reversal behavior of Pq. Note that there exist orbits connecting P 4 and P 5 with Pq 
and P 7 , which implies that Q 0 can cross zero, i.e. H = 0, during the evolution, and thus 
the universe exhibits a bounce or a cosmological turnaround. Finally, the system admits an 
static solution, Pg, which always behaves as a saddle point, hence it cannot represent the 
late-time universe. We mention that the above features are not obtained in the flat or open 
curvature, where H cannot change sign. 

In order to present the above behavior in a more transparent way, we evolve numerically 
the cosmological equations and in Fig. 2 we depict the corresponding phase-space behavior. 


- 19 - 







1 


0 


DE 


0.8 

1 V 

\ 










s’ 

/// 


1 \\ 

A 

A 


\ 


.. ~ 

~ 



s* 

s' 

/ 


V 

/ 


i \ 

V 

V 

\ 

O 


■.... 




S’ 

A, 

'/A 


/ 



1 ■ 





v > 

. 




S/ 



/ 


t 

0.6 

1 . 

l 



A 








ft 

/ 


t 

i 

1 

!\ 


\ 


\, ^ 



s* 

// 


t 

/ 

/ 

t 


1 

i 

\ 


\ 

V 

v\ 



Sj 

// 

'? 

t 

/ 

/ 

t 


1 

l 


i 








t 

t 

f 

t 

t 

fl A- 

1 

1 

1 

i 

l 


Iv 



7 

/ 

t 

t 

f 

t 

t 

U.H- 

1 

i 

1 


1 

1 

i) l 



f 

t 

f 

t 

1 

t 

t 


1 

1 

1 

i 

l 

l 

IN 



\ 

[ 

t 

t 

t 

t 

t 


1 

l 

1 

H 

1 

l 

/ / 





\ 

\ 


t 

I 

0.2 

1 

l 

/ 

1 

/ 

A 




V* 

V 

A 

\ 

\ 

\ 

t 


i 

1 

/ 

/ 

/. 






88 


A 

\ 

\ 





0.5 


Figure 2. The phase-space behavior of time asymmetric cosmology of Model I: f(V) = g\V m , with 
positive curvature, p = 0.9 (i.e. m = —0.033,), and w m = 0. In this specific example the universe 
is led to either (a) the dark-energy dominated, accelerating solution P$ or (b) the matter dominated, 
contracting solution P§. 


As we can see, in this specific example the universe results in the dark-energy dominated, 
accelerating solution P 5 or in the matter dominated, contracting solution Pq. 

4.2 Model II: f(V) = g 2 e xv 

In the case where f(V) = g 2 e xv , i.e. when g(a) = g 2 ^e Xa3 , with g 2 a constant and A 
a parameter, with open or zero curvature, and A > 0, the scenario at hand exhibits five 
isolated critical points and one curve of critical points presented in Table 3. Amongst them, 
only point Qq behaves like a stable one (although nonhyperbolic) and thus it can be the 
late-time state of the universe. It corresponds to a dark-energy dominated universe, in which 
the dark-energy equation-of-state parameter lies in the phantom regime. Note however that 
as the universe approaches this point, the deceleration parameter q decreases monotonically, 
resulting to a divergence at the critical point. In particular, as the scale factor increases 
and the dark energy term becomes dominant, we can obtain an approximate solution for the 

scale factor, namely the inverse of t — to = E ^ 3g “, ^ = e~ a3x O where to = — c\/g 2 , 

with Ei (z),z < 0, the exponential integral function and c\ an integration constant, and we 
can immediately see that the scale factor diverges at a finite time, which is the realization 
of a Big Rip [ 68 ]. This behavior was expected, since for A > 0 the extra, time-asymmetric, 
term that constitutes the effective dark energy sector increases monotonically. Hence, for 
these parameter choices, the scenario at hand does not correspond to the usual classes of 
cosmological models, and thus it should not be considered as a successful one. In Fig. 3 we 
depict the phase-space behavior of such a scenario, arising from numerical elaboration. As 
we observe, in this example the universe results in the dark-energy dominated, accelerating 
solution Qq. 
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Figure 3. The phase-space behavior of time asymmetric cosmology of Model II: f(V) = g 2 e xv , with 
negative curvature, w m = 0 and A > 0 (the specific value of A is not relevant, only its sign, since 
it has been absorbed into the auxiliary variable T according to (3.21)). In this specific example the 
universe is led to the dark-energy dominated, accelerating solution Qq. The bold dashed line named 
Q 5 in general presents saddle behavior, however it is a local source for all the orbits (which have the 
shape of straight lines connecting it with Qq) located at the invariant set T = 1. 


In the case of zero or open curvature and A < 0, the model exhibits five isolated critical 
points and one curve of critical points, displayed in Table 4. Amongst them, point Qio 
behaves as stable for the flat models, and thus it can attract the universe at late times. 
However, it corresponds to a dark-matter dominated universe, and thus it is not favored 
by observations. This was expected, since for A < 0 the effective dark-energy terms are 
redshifted away in a much faster way (due to the exponential) than the matter contribution, 
leaving the universe matter dominated. Nevertheless, one could improve this behavior by 
the addition of an explicit cosmological constant, in which case he could get the correct 
thermal history, namely the succession of matter and dark-energy eras. However, since in 
this work we are interested in investigating the effects of the pure time-asymmetric cosmology, 
without the explicit presence of a cosmological constant, we do not examine such a possibility 
further. Additionally, as we describe in detail in Appendix A.2, the nonhyperbolic curve of 
critical points Qn, with the exception of its endpoint with = 1, behaves as saddle. For 
0 < Jlfc < 1, it corresponds to a universe with Qde = 0, however not completely matter- 
dominated, since the curvature contribution remains non-zero. Another interesting point 
located on the curve Qu is the one corresponding to complete curvature domination, namely 
with = 1. This point is indeed a stable late-time state of the universe. Similarly to 

Qio, the above features are not favored by observations to be the late-time state of the 
universe, however these curves of points could be a good candidate for the description of its 
intermediate phases, especially under the addition of an explicit cosmological constant. In 
Fig. 4, through a numerical elaboration, we present the phase-space behavior of this model. 
As we see, in this example if the universe starts with ttf. = 0 it results in the dark-matter 
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dominated solution Qiq. On the other hand, if 0*. > 0 initially then the universe results in 
the curvature-dominated solution (fi*,, Qde, T±) = (1,0,1) located on the bold dashed line 
Qn- 



Figure 4. The phase-space behavior of time asymmetric cosmology of Model II: f(V) = g 2 e xv , with 
negative curvature, w m = 0 and X < 0 (the specific value of X is not relevant, only its sign, since 
it has been absorbed into the auxiliary variable Tj according to (3.31)). In this specific example the 
universe is led to the dark-matter dominated solution Q io (if life = 0 at the initial state), or to the 
curvature-dominated solution located on one endpoint of line Qn, namely (Xik,XlDE,Ti) = (1,0,1) (if 
Clk > 0 at the initial state). All other points of the curve Qn, which is represented by a bold dashed 
line, behave as saddle. 

In the case of positive curvature and A > 0 the scenario at hand exhibits four phys¬ 
ical critical points, and one curve of critical points, namely Qie, which is the straight line 
joining the points (0,0,1) and (1,0,1) (with the left endpoint not included), corresponding 
to expanding cosmologies, which are summarized in Table 5. All these critical points have 
contracting partners via the discrete symmetry (3.42), which are displayed in Table 6. Addi¬ 
tionally, there exists a line of static solutions namely Si, however since they are saddle they 
cannot attract the universe at late times. Amongst all these points, the late-time attractors 
are the expanding solution Q 17 and the contracting R13. In particular, Qn corresponds to 
a dark-energy dominated universe in which the dark-energy equation-of-state parameter is 
phantom-like. Note however that as the universe approaches this point, the deceleration 
parameter q decreases monotonically, resulting to a divergence at the critical point. Using 
similar arguments as for point Qq for the open or zero curvature case, it can be shown that 
it is of a finite-time type, namely a Big Rip [68]. Similarly to the open or zero curvature 
case, this behavior was expected, since for A > 0 the extra, time-asymmetric, term that 
constitutes the effective dark energy sector increases monotonically. Additionally, there is 
another stable late-time solution, namely the matter-dominated point A43 which ends in a 
Big-Cruch. In Fig. 5 we depict the phase-space behavior of this scenario. As we see, in this 
example the universe results in the dark-energy dominated, accelerating solution Q 17 or in 
the Big-Crunch singularity A43. 
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Figure 5. The phase-space behavior of time asymmetric cosmology of Model II: f(V) = g 2 e xv , with 
positive curvature, w m = 0 and A > 0 (the specific value of X is not relevant, only its sign, since it has 
been absorbed into the auxiliary variable T according to (3.27)). In this specific example the universe 
is led to either the dark-energy dominated, accelerating solution Qn, or to the contracting solution 
i?i3. The dot-dashed (red) line represents the curve of static solutions Si. Notice the presence of orbits 
crossing the line Q = 0, i.e. H = 0, which correspond to transitions from expanding to contracting 
cosmologies and vice versa, that is to cosmological turnarounds and bounces. 


In the case of positive curvature and A < 0, the model exhibits four isolated critical 
points and one curve of critical points, corresponding to expanding cosmologies, displayed in 
Table 7. Each of the above critical points have contracting partners via the discrete symmetry 
(3.42), which are displayed in Table 8. Amongst them, point behaves as stable, and thus 
it can be the late-tinre state of the universe. However, it corresponds to a contracting dark- 
matter dominated universe, and therefore it is not favored by observations. Similarly to 
the open or flat case, this was expected since for A < 0 the effective dark-energy terms are 
redshifted away in a much faster way than the matter contribution. Additionally, there exists 
a line of static solutions which are saddle, while, as we describe in detail in Appendix A.4, the 
nonhyperbolic curves of critical points Q 20 and Q 21 behave typically as saddle. In Fig. 6 we 
present the phase-space behavior for the model at hand, where we observe that the late-time 
attractor is the contracting solution R\g. Additionally, the figure shows orbits exhibiting the 
crossing of the Qo = 0 line, which correspond the transition from contracting to expanding 
cosmologies and vice versa. 

5 Conclusions 

In this work we studied the cosmological behavior in a universe governed by time asymmetric 
extensions of general relativity. This novel modified gravity is based on the addition on 
the Hamiltonian framework of new, time-asymmetric, terms, in a way that the algebra of 
constraints and local physics remain unchanged [38]. However, at cosmological scales these 
new terms can have significant effects that can alter the universe evolution, both at early and 
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Figure 6. The phase-space behavior of time asymmetric cosmology of Model II: f(V) = g 2 e xv , with 
positive curvature, w m = 0 and X < 0 (the specific value of X is not relevant, only its sign, since it has 
been absorbed into the auxiliary variable T\ according to (3.31)). In this specific example the universe 
is led to the contracting solution Ris- Additionally, the figure shows orbits exhibiting the crossing of 
the Q o = 0 line, which correspond the transition from contracting to expanding cosmologies and vice 
versa. 


late times. In particular, assuming that the new terms in the Hamiltonian are proportional to 
an arbitrary function of the spatial volume, we finally obtain modifications of the Friedmann 
equations depending on an arbitrary function of the scale factor. Definitely, the capabilities 
of such cosmological constructions are huge. 

We considered two basic ansatzes for the aforementioned modification, namely a power 
law and an exponential one. We mention that we did not consider an explicit cosmological 
constant, since we desired to investigate the pure effects of the new terms. In order to bypass 
the complexity of the equations, we applied the dynamical systems method, which allows 
to reveal the global behavior of time asymmetric cosmology, independently of the details 
of the evolution and the specific initial conditions. In particular, we extracted the critical 
points of the scenario and we examined which of them are stable and thus they can be the 
late-time state of the universe, calculating also the corresponding observables, such as the 
various density parameters and the deceleration parameter. 

For the power-law ansatz we found that the universe can result in a dark-energy domi¬ 
nated, accelerating universe, where the dark-energy equation-of-state parameter wde can lie 
either in the quintessence or in the phantom regime, or even behave as an effective cosmo¬ 
logical constant giving rise to a de Sitter universe. Moreover, by suitably choosing the model 
parameter, one can obtain a wp>E in agreement with observations. 

For the exponential ansatz we showed that for positive exponential coefficient at late 
times the universe is attracted by a dark-energy dominated universe, in which wde lies in 
the phantom regime, resulting finally to a finite-time Big-Rip singularity (due to the expo- 
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nential increase of the novel terms). On the other hand, for negative exponential coefficient 
the universe results to a dark-matter dominated universe (due to the exponential decrease 
of the novel terms comparing to the matter sector), which is not favored by observations. 
Nevertheless, one could improve this behavior by the addition of an explicit cosmological 
constant, in which case he could get the correct thermal history, namely the succession of 
matter and dark-energy eras. Finally, note that in the case of closed curvature, the universe 
may experience a cosmological bounce or turnaround, or even cyclic behavior. 

Concerning phenomenology, we should mention that in the scenario at hand the left 
handed neutrinos propagate differently than the photons [38], since the latter propagate 
according to the usual connection of the spacetime metric, while the former propagate ac¬ 
cording to the Ashtekar connection and geometry. Hence, if one desires to be in agreement 
with observations, for instance with the data from SN1987A supernova which show that 
massless neutrinos propagate similarly to photons with an error less than 1CP 9 [69-71], then 
he should impose the new time-asymmetric modifications to be small, as expected. Interest¬ 
ingly enough, even if one considers the extreme realization of the above requirement, namely 
to assume that the new terms tend asymptotically to zero (instead of being increasing) as the 
universe expands, one can still have significant effects at large scales, that can radically alter 
the universe behavior (for instance in the power-law modification with f(V) = giV m and 
p = 3m + 1, for the parameter window 0>m>—1/3 one has an asymptotically vanishing 
modification term which is nevertheless able to drive late-time acceleration (since 1 > p > 0) 
since its tends to zero as a( 2p ~ 2 ' 1 i.e slower than the matter and curvature contributions in the 
Friedmann equation). Hence, one can easily pass all the cosmological tests, and definitely 
all the Solar System ones. An interesting study would be to examine the bounce realization, 
since in such a case one would expect the time asymmetry to lead to distinguishable signa¬ 
tures on observations, especially having in mind the different behavior of the spacetime and 
Ashtekar related quantities. Additionally, an important and necessary investigation would 
be to examine the cosmological perturbations and their relation to various observables, either 
at early, inflationary times, or at late epochs. Since both these studies lie beyond the scope 
of the present work they are left for future projects. 

In summary, the cosmological application of time asymmetric extensions of general 
relativity has many capabilities and thus it can be a good candidate for the description of 
the universe, that is worthy to be studied further. 
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A Stability of the nonhyperbolic critical points of Model II: f{V ) = g 2 e AV 

In this Appendix we investigate the stability of the nonhyperbolic critical points that appear 
in the analysis of Model II in subsection 3.2, using the center manifold method [61], since in 
this case the simple linear analysis is not adequate. 
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A.l Zero or negative curvature and A > 0 

In the case of zero or negative curvature and A > 0, we extract two isolated nonhyperbolic 
critical points, and a curve of nonhyperbolic critical points, displayed in Table 3. Since point 
Q 4 and the curve Q§ have at least one unstable eigen-direction they will definitely be non¬ 
stable (i.e. saddle or unstable), and hence we do not need to perform the center manifold 
analysis, since in this work we are interested in the stable late-time solutions. Thus, we 
restrict our analysis in the case of Q&. 

We introduce the new variables 


e = l-T, x = Q k , y = l-^DE, (A.l) 

in order to translate Qq to the origin, and thus we obtain the system 

| = 3{f - l)f 2 , (A.2a) 

dx 

— = -x {e(3 w m x + x - 4) - 3y [(w m - l)e + 2] + 6 } , (A.2b) 

dr] 

^| = (1 -y) {(3 w m + l)xe - 3 y [(w m - l)e + 2 ]} , (A.2c) 

where the local center manifold of the origin (e, x, y ) = ( 0 , 0 , 0 ) is tangent to the e-axis. 
Hence, it can be written locally as the graph 

{{e,x,y) : x = h x (e),y = h 2 (e),hi(0) = 0,/i 2 (0) = 0, M°) = 0,/i' 2 (0) = 0, |e| < T}, (A.3) 


where 5 is a suitably small number. The functions hi and h 2 must satisfy the quasilinear 
system of differential equations 

3(e - 1 )e 2 h'i(e) + h x {e) {e [ 3w m hi(e ) + h x (e) - 4] - 3 h 2 (e) [(w m - l)e + 2] + 6 } = 0, (A.4a) 
[1 - h 2 (e)]hi(e)(3w m + l)e + 3 h 2 (e) [( w m - l)e + 2)] - 3(e - l)e 2 ti 2 (e) = 0. (A.4b) 

This system admits the following solutions: 

1 . the point: 



hi(e) = 0 , 


(A.5) 


h 2 {e) = 0 , 


(A. 6 ) 

2 . the 1 -parameter solution: 



hi(e) = 0 , 


(A.7) 

f— , -- , e /0 

h 2 (e) = < e cl+ « (l-e)“" l + 1 +e' 1 " m+1 , 

l 0 , e = 0 


(A. 8 ) 

3. the 2-parameter solution: 




e c 2 £ 2/3(l e ) u '>T l + ^ 

e ^0 

e = 0 


= < 

1 e c 2e 2 / 3 (l—e)^ m + ^+cie 2 / € (l— e) Wrn + 1 -\-e w ™'+ 1 5 

l 0 , 

(A.9) 

1 

e =2 e 2 /3(l e ) w »w + 3| e U'm+l 

e^° 

e = 0 


h 2 (e) = i 

e c 2e 2 / 3 (l—e)^ m "*”^+cie 2 / € (l— e) Wrn + 1 -\-e w 'm+ 1 ? 

l o, 

(A. 10 ) 
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These three classes of solutions satisfy the smoothness conditions required in order to obtain 
the center manifold of the origin (note that the expression for the center manifold is not 
unique). Thus, we conclude that the evolution on the center manifold is given by the equation 


which admits the solution 


de 

dfj 


-3(1 -e)e 2 , 


>7 = c 1 + i 


- + 2 tanh 1 (1 — 2 e) 


Cl + 3^ 


log(e) 

3 


e 

3 


+ 0(e 2 ), 


(A.ll) 


(A-12) 


and therefore by inverting the above expression we find e(fj). It is easy to see that e -> 0 as 
77 —» 00 and that e —» 1 as fj —>• — 00 . Hence, we deduce that the center manifold of Qq is 
stable [61]. 


A.2 Zero or negative curvature and A < 0 

In the case of zero or negative curvature and A < 0, we extract two isolated nonhyperbolic 
critical points, and a curve of nonhyperbolic critical points, which are presented in Table 4. 
Since point Q 12 has at least two unstable eigen-directions it will definitely be non-stable, and 
hence we do not investigate it further. 

In order to examine the stability of Q 10 using the center manifold theorem we introduce 
the variables 


e = 1 — Ti, u = Qde, v = Qk, 

(A.13) 

with evolution equations given by 


de q ( -1 \ 2 

dij ~ 3(6 " 1)6 ’ 

(A. 14a) 

d'U r r 

— = — u {e 3 u> m (u + v — 1) + 9u + v — 6u — 9e + 6} , 
dij 

(A.14b) 

dv 

-TZ = -V {3 u [(w m + 3)e - 2] + (v - 1)(3 w m + l)e} . 
dij 

(A.14c) 

The center subspace of the origin of (A.14) is spanned by the vectors (1,0,0) and (0,1,0), 
which implies that the local center manifold of the origin can be written locally as the graph 
{(e, u, v) : v = h(e,u),h( 0,0) = 0, Dh(0,0) = 0 , |(e,u) < <5}, where Dh is the matrix 

of derivatives, 5 is a suitably small constant, and h(e, u ) satisfies the quasilinear partial 
differential equation 

r)h r)h 

u— {(3w m +)eh + 3 (u - 1) [(w m + 3)e - 2]} - 3(e - l)e 2 ^- 
ou oe 


- h{(3w m + 1 )e{h — 1) + 3u [(w m + 3)e — 2]} = 0. 

(A.15) 


Assuming that h(e,u) = uf(e) and lim e _>.o/(e) = lim^o f'( e ) = 0, and substituting in 
(A. 15), we obtain 


3(e - 1 )e 2 /'(e) + ( 8 e - 6 )/(e) = 0 , 


(A.16) 
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which has the general solution 


f(e) = 


c±e 


2A,2/3 


(A.17) 


(1 — e) 2 / 3 5 

and the trivial solution f(u ) = 0. However, the general solution leads to lim e ->o/(e) — 
sgn(ci)oo, lim^o/(e) = —sgn(ci)oo, and hence it does not satisfy the imposed limits. Thus, 
the only accepted solution is the trivial one, which implies h(e, u ) = 0. Hence, for this case 
the dynamics on the center manifold is governed by 

= 3(e — l)e 2 , (A.18a) 


df) 

du 

df) 


= -3 (u - 1 )u \{w m + 3)e - 2] 


Eliminating time and integrating out we finally acquire 

1 


u{e) = 


g c l + 7 e «J m + l(q _ g'j— Wm—l _|_ \ 


(A. 18b) 


(A.19) 


which satisfies u —>• 0 as e —>• 0. This feature implies that Q io attracts the orbits contained in 
its center manifold (that is the 2D set T-Q.ee), and thus this nonhyperbolic point behaves 
as stable. 

In order to examine the stability of the curve of critical points Qn (with Qkc € (0,1]), 
using the center manifold theorem, we introduce the variables 


e = 1 — Ti, u = f2fc c (l — Qde) ~ v = ^de, 
which satisfy the evolution equations 

^ = -3(1-^. 
df) 


du 

df) 


dv 

di) 


= e {u 2 (3w m + 1) + u [-3 v(w m + 3) + 3w m + 1]} 

+ H fcc e{tt [(u - 2)(3 w m + 1)] + (1 - H fcc )(u - l)(3w m + 1)} + 6 uv, 
= e [uv{3wm + 1) + 3v(l - v)(w m + 3)] 


(A.20) 


(A.21a) 


(A.21b) 


(A.21c) 


+ H fcc e [v(v - 1)(3 w m + 1)] + 6v(v - 1). 

Since the center subspace of the origin of (A.21) is spanned by the vectors (1,0,0) and 
(0,1,0), we deduce that the local center manifold of the origin can be written locally as the 
graph {(e, u,v) : v = h(e, u ), h(0, 0) = 0, Dh(0, 0) = 0 , ||(e, u)\\ < <5}, with 6 a suitably small 
constant, and where the function h(e, u) that defines the center manifold must satisfy the 
quasilinear partial differential equation 

h{(3w m + 1 )(Hfcc - l)^fc c e - u{e [3 w m (Q k c - 1) + ^kc - 9] + 6} | 


e(3 Wyyi -{- 1)(m Hfc c )(w H kc T I) 


dh 

— 3(e — l)e 2 ——b {e [3u> m (H/ cc — 1) + — 9] + 6} h 2 

oe 

+ h {e [(it — Hfc c )( 3w m + 1) + 3(w m + 3)] — 6} = 0. 


(A.22) 
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For Q kc / l,w m / —1/3 the above equation should be integrated numerically. 

We will proceed using Taylor expansion. In particular, the solution /i(e, u) must satisfy 
the conditions h(0, 0) = 0, Dh(0,0) = 0 , that is it must be at least of second order in the 
variables e and u. Hence, we assume that h(e,u) = one 2 + a^ew + a, 2 2 u 2 + 0( 3), where 
0(3) denotes terms of third order on the vector norm, i.e. terms like e 2 u, eu 2 , e 3 , u 3 . These 
terms and higher-order terms neglected in the approximation scheme. Substituting back this 
expression for h, neglecting third-order terms, comparing terms of the same power, equating 
to zero the coefficients, and assuming that Q k c / 1 / —1/3, we obtain that a good 

approximation of the center manifold is given by 


h(e, u) = ^ai 2 (l + 3n; m )(l - H fcc )H fcc e 2 + ai 2 eu + 


3012 


f ^ \ - ' >i*n- -ivw-™- ' '(-1,0 n U ' (A.23) 

6 (1 T 3w m )y\ ^kcj^kc 

Therefore, we deduce that the dynamics on the center manifold is determined up to third 
order by 

de 
df] 

+ l)e [(Q kc - l)fi fcc - u{ 2n kc - l)] . (A.24b) 


= —3e 2 


(A. 24a) 


du a 

-JZ = \3w% 

di1 


The system (A.24) admits the general solution 

1 


e(v) = o- 

3 rj - ci 

u(fj) = c 2 (3 f, - Cl ) -i(3» ra +i)(2n hc -i) + 

Z\lkc — 1 


(A.25) 
(A.26) 


Observe that as r) —> +00, e — > 0, but u departs from zero and becomes unbounded in the case 
^kc < \i or tends to ^20/ for tt kc > 1 as f —>• +00, which is nonzero since f l kc ^ {0,1}. 
Thus, the origin is unstable along the u-axis and stable along the e-axis. Summarizing, the 
line of fixed points Qu behaves as saddle, provided that Q kc 7^ l,w m / —1/3. 

Let us mention that the above analysis is essentially an approximation. Nevertheless, 
there is a special point of the curve Qn, namely (fifc, Qde, 7i) = (1)0,1), for which the 
above procedure is not valid, that allows for an analytical application of the center manifold 
analysis. It corresponds to Q kc = 1 in (A.21). Setting £l kc = 1 in (A.22) we obtain the 
simpler quasilinear partial differential equation 

r)h r)h 

- u [(6 - 8 e)h + (u- 1)(3 w m + l)e] — - 3(e - l)e 2 — 

+ h [(6 — 8 e)h + e(3 uw m + u + 8) — 6] = 0. (A.27) 

Given the solution v = h(e,u), the dynamics on the center manifold is determined by 


de 

df 

du 

df 


= -3(1 - 6)6 2 

= u[(u — l)(3w m + l)e + h(e, u)( 6 - 8e)] . 


(A.28a) 
(A.28b) 


Assuming that h(e,u ) = uf(e) and lim e _>.o/(e) = lim e _»o f'(e) = 0, and substituting into 
(A.27), we obtain 


{/(e) [{w m + 3)e - 2] - (e - l)e 2 /'(e)} = 0, 


(A.29) 
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which has the general solution 


f(e) 


cie" 2 / e (l - e ) Wm+ 1 e~ Wrn ~ 1 , 

0 , 


e/0 

e = 0 ’ 


(A.30) 


which indeed satisfies the imposed limits. Hence, the dynamics on the center manifold is 
governed by the evolution equations 


de 

df 


-3(1 - e)e 2 



c lU e ~ 2/e (6 - 8e)(l - e )™-+ 1 e -“'"‘- 1 + („ 


1)(3 w m + l)e 


(A.31a) 
(A.31b) 


Eliminating the time variable we find that the system (A.31) can be expressed as 


with //(e) 


3(1 - e)e 


du(e ) 
de 


2cie -2 / f (4e — 3)(1 — e ) Wm+1 


= (2w m + l)u(e) + //(e)«(e) 2 (A.32) 

e -w m -2 _ 2 >w m — 1, which admits the quadrature 


u{e) 


(1 -e)- 




c 2 -/^(e)(l-e) 


w m Ig'U'm § 


(A.33) 


Since //(e) —>• — 3w m — 1 as e —>• 0, we can integrate the above quadrature in the approximation 
e —> 0, obtaining 


u(e) 


1 

c 2 (l - e ) Wm+ ^+ 1 


(A.34) 


which tends to zero as e —>• 0, for w m > —4- Hence, we deduce that the center manifold 
associated to the point (1, 0,1) is stable. Indeed, this behavior is the typical one for w m > — 
as can be verified by Fig. 4. 


A.3 Positive curvature and A > 0 

In the case of positive curvature and A > 0, we extract two isolated nonhyperbolic critical 
points, and a curve of nonhyperbolic critical points corresponding to expansion, which are 
presented in Table 5. Each of the above points/curve in Table 5 has a partner through the 
symmetry (3.42), which represents a contracting cosmology, and are displayed in Table 6. 
Amongst them in this Appendix we analyze only the nonhyperbolic fixed points that might 
be late-time attractors (for instance points like Q 15 and the curve of critical points Q that 
have at least one unstable eigen-direction will definitely be either unstable or saddle and 
thus we do not investigate them further). These are the expanding solution Qn and the 
contracting solutions R 15 and R\q. We remind that the points Qi and their contracting 
partners points Ri through the symmetry (3.42), exhibit opposite dynamical behaviors, and 
thus from the following analysis we also obtain information for the contracting solution Rn 
and the expanding ones Q15 and Q±q. 

In order to examine the stability of the contracting solution A45 we introduce the vari¬ 
ables 


e = l-T, x = l + Q 0 , V = @de, 


(A.35) 
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and therefore the autonomous system (3.40) becomes 


de 

df 

dx 

df 

dy 

df 


3(1 - x)(e - l)e 2 , 

(A.36a) 

i(x - 2)x(e(3 w m (y - 1) - 3y - 1) + 6 y), 

(A.36b) 

3(1 - x)(l - y)y[(w m - l)e + 2], 

(A.36c) 


The center manifold of the origin of (A.36a) is spanned by the vectors (1,0,0) and (0,1,0), 
which implies that the local center manifold of the origin can be written locally as the graph 
{(e,x,y) : y = h(e, x), h(0, 0) = 0, Dh(0,0) = 0, ||(e, u)|| < <5}, with 5 a suitably small 
constant and Dh the matrix of derivatives. The function h(e, x) satisfies the quasilinear 
partial differential equation 



dh r 

2)x— 3(w m - l)e + 2 h — (3w m + l)e - 3(x 
ox L J 

— 3(x - l)[(tc m - l)e + 2 ](h — 1 )h = 0. 


D(e-1)6 : 


,dh 

Te 


(A.37) 


The equation (A.37) admits the solutions: 
1. The trivial solution h(e,x) = 0, 


e 2 / £ (e—l)(l—e) Wm 


2. the one-parameter solution h(e,x) = { -e c ie w m+^- e 2 / £ {\-e) w m-\-e 2 / e e{i-€) w m e ^ ® . 

[ 1 e = 0 

Only the trivial solution satisfies the conditions h( 0, 0) = 0, Dh(0, 0) = 0. Henceforth, the 
dynamics on the center manifold is governed by 


= 3(1 — x)(l — e)e 2 , 
dx l, 

— = -(3 w m + l)(ic - 2 )xe. 
ar z 


(A.38a) 
(A.38b) 


Eliminating the time variable, r, and using the chain rule for derivatives we find that the 
orbits on the invariant manifold satisfy 


x'(e) = 


(3w m + l)(x(e) - 2)x(e) 
6(e - l)e(x(c) - 1) ’ 


which admits the general solutions 


r ^j e Wrn+ ^ — e 2ci (1 -e)^ m+ ^. 


(A.39) 


(A.40) 


None of these solutions satisfy the condition x(0) = 0, indeed x becomes infinity as e —>• 0. 
Thus, any solution starting with e / 0 and x ^ 0 departs from the origin along the x-direction, 
which implies that R\^ behaves as a saddle. 

In order to examine the stability of the contracting solution we introduce the vari¬ 
ables 

2 1^0 1 

ui = Qo --—-1- Qoc, U 2 = - (Qoc 2 — l) (1 — T)(3w m + 1), v = Qde, (A.41) 

ZU() C z 
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where Q q c is a constant Q q c € (0,1), which satisfy the equations 


du\ 

df 

1 

~~ 2 

dU 2 

df 

dv 

df 


3{v - l)v [ Q 0c 2 (v - 2) + 2Q 0c m - v] [w m (3 Q 0c 2 + u 2 - 3) + Q 0 c 2 - u 2 - l] 


2Qoc (3w m + 1) 


QocV , V 

~ + 2 Qf + Qac ~ l] ~ 


2u2(3v(w m - 1 ) - 3 w m - 1 ) 
(Qoc 2 - l) (3w m + 1 ) 


+ 6v 


3 u 2 2 [Qoc 2 {y - 2) + 2Q 0c ui - v ] [Qoc 2 (-(3w m + 1)) + 2 u 2 + 3u; m + l] 


Qoc (Qoc 2 -lj (3w m + iy 


3{v — l)n 


Q 0c v 


2(5oc 


+ Qoc — 


2 u 2 (w m - 1) 


(Qoc~ — l) (3 w m + 1) 


+ 2 


(A.42a) 
(A.42b) 
(A.42c) 


Note that by definition u 2 < 0. 

The center subspace of the origin of (A.42) is spanned by the vectors (1,0,0) and 
(0,1,0), which implies that the local center manifold of the origin can be written locally 
as the graph {(ui,rt2,u) : v = h{e,u),h{ 0,0) = 0, Dh(0,0) = 0 , ||(iti, u 2 )\\ < 5}, with 5 a 
suitably small constant and Dh the matrix of derivatives. The function h{ui,u 2 ) satisfies 
the quasilinear partial differential equation 


3 u 2 2 [3 (Qo 2 - 1 )w m + Qoc 2 - 2 u 2 - l] [(Qoc 2 -l)h + 2Q 0c (u 1 - Q 0c )] dh 


Qoc {Qo 2 -l) {3w m + iy 


du 2 


1 dh 

2 du\ 


{Qo 2 ~ 1 )h 
2Qq c 


Qoc T 


-1 > [ 2U2 \ 3(w ™ 3Wm 1] +6 h 

\ {Qoc - l) (3 Wm + 1) 


3 [w m (3Qo c 2 + U 2 - 3) + Qo 2 -U 2 - l] [h - 1 )h [{Qo 2 -l)h + 2Q 0c {ui - Q 0c )] 


Qoc (3w m + 1) 


+ 3 


2 U 2 {w m - 1) 


(Qoc 2 - 1) {3w m + 1) 


+ 2 


{h — 1 )h 


{Qo 2 ~ 1 )h 
2 Qqc 


+ Qoc — Ul 


= 0. 


(A.43) 


Assuming that h is locally given by h = anuf + ai 2 uiu 2 + 022^2 + 0(3), where 0(3) denotes 
terms of third order in the vector norm, plugging back into (A.43), comparing equal powers 
in the variables u\ and u 2 and equating to zero the corresponding coefficients, we obtain 
up to third order that an = —3ai2<5oc and a 22 = — That is, the graph of the center 
manifold is given up to third order by v = — 3ai2Qoc^f + 012^1^2 — -£^ u 2 - Therefore, by 
neglecting the third-order terms we find that the dynamics on the center manifold is governed 
by equations 


du\ 

df 

du 2 

df 


u 2 + 


2Qoc^1^2 

1 - Qgc ’ 


^Qoc ^2 

(l — Qq c ) (3 w m + 1) 


(A.44a) 
(A.44b) 
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Integrating (A.44) it follows 


ui(t) 


U 2 (t) 


2c 2 Q 0c [ci ( Ql c - l) (3w m + 1) + 6 Q 0 c t] ^ m+ s + Q\ c - 1 

2 Qqc 

_ (Qoc ~ -Q ( 3w m + 1)_ 

ci (Qlc ~ !) (3w m + 1) + 6Q 0c t 


(A.45a) 
(A.45b) 


Taking the limit as r —>• oo in the above expressions we obtain {u\. u 2 ) —>• (c 2 oo,0),c 2 ^ 
0, U 2 / 0. In the special case c 2 = 0 we obtain the limits (rq, tt 2 ) —» ( 2 q 0 >0)- In both cases 
the origin is unstable along the tq-axis. Since it is stable along the tq-axis, it follows that 
Riq is a saddle. 

Finally, let us examine the stability of Qn using the center manifold theorem. We 
introduce the variables 


e = 1 -T, x = 1 - Qq, y = 1 - Q DE , (A.46) 

and therefore the autonomous system (3.40) is equivalent to the system 

— = —3(1 — x)(l — e)e 2 , (A.47a) 

dr 

dx 1 

— = -~(x - 2)x [3 y{{w m - l)e + 2} + 4g - 6] , (A.47b) 

ar Z 

= -3(1 - x)(l - y)y \(w m - l)e + 2] . (A.47c) 

The local center manifold of the origin (e, x, y) = (0,0, 0) is tangent to the e-axis. Thus, 
it can be written locally as the graph {(e,x,y) : x = /q(e),y = /i 2 (e), /q(0) = 0,/i 2 (0) = 
0 , h\ ( 0 ) = 0 , ^ 2 (0) = 0, |e| < 5}, with 5 a suitably small number. The functions /q and h 2 
must satisfy the quasilinear system of differential equations 


[Me) “ 2]/q(e)[3/i 2 (e)((u; m - l)e + 2) + 4e - 6] - 6(e - l)e 2 (/q(e) - 1 )h[(e) = 0, (A.48a) 
M e ) - 1] [(/i 2 (e) - 1 )h 2 (e)((w m - l)e + 2) - (e - l)e 2 h' 2 (e )] = 0, (A.48b) 

which admits the general solution satisfying the conditions /q(0) = 0, /i 2 (0) = 0,^(0) = 
0,14(0) = 0, namely: 


/11(e) 


/12(e) 



\f~ 


e 2c 2e 2/3(l e )“m + ^| e ci + | + + 1 


V eci+ ^(i — f^j w m +1 -J-g'ie'm +1 

0, 


e/0 


e = 0 




g 0 ! - ! - e —e^m+l-j-e^rn+l 

0, 


e/0 
e = 0 ’ 


(A.49a) 

(A.49b) 


where ci and c 2 are integration constants, as well as the trivial solution /q(e) = 0, /i 2 (e) = 0. 
Note that the expression for the center manifold of the origin is not necessarily unique. 

For the expression of the center manifold of the origin given by (A.49), the dynamics 
on it is given by 


de 

df 


-3/(e)(l-e)e 2 , 


(A.50) 
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where 


f(e) = 


g2c 2e 2/3(l _ e )Wm+5 g c l + 7 (1 — fjHm + 1 + e Wm + l 


^/e ci+ ?(l - e)™- +1 + e“ m +1 
Since /(e) > 0, the flow of (A.50) is equivalent to the flow of 

de 


dt, 


= -3(1 - e^ 


where we have introduced a time rescaling. The general solution of (A.52) reads 


£ = c 3 + 


1 


1 


+ 2 tanh *(1 — 2e) 


= C 3 + 1 _ Ml) _ i + o (e 2 ) . 
3e 3 3 y 1 


(A.51) 


(A.52) 


(A.53) 


Since e —>• 0 as £ —>• oo, the center manifold of Q17 is stable, and it corresponds to the 
late-time attractor. Additionally e —>• 1 as £ —>• —00. Note that the relation with the original 
time variable is obtained through the quadrature 


r = 


?(e)de 1 
f(e) 


\J e Cl+ £ (1 - e)"' m+1 + e w m +1 


: de. 


3 d (e - l)e 2 -y/e 2c2 e 2 / 3 (l - e) Wm+ i + e Cl+ ? (1 - e)' u '" i+1 + e"'" i+1 

(A.54) 

If the center manifold is given by the trivial solution h\{e) = 0,/12(e) = 0 we deduce 
that the evolution on it is dictated by 


de -on \2 
---3(l-e)e , 


which admits the solution 


1 

f = C1 + 3 


—f 2 tanh 1 (1 —2e) 




(A.55) 


(A.56) 


Since e —>• 0 as r —>• 00, the center manifold of Qn is stable, and it corresponds to the 
late-time attractor. 


A.4 Positive curvature and A < 0 

In the case of positive curvature and A < 0, we extract two isolated nonhyperbolic criti¬ 
cal points, and a curve of nonhyperbolic critical points, representing expanding solutions, 
which are presented in Table 7. Each of the above points/curve in Table 7 has a partner 
through the symmetry (3.50), which represents a contracting cosmology, and are displayed 
in Table 8. Amongst them we analyze only the nonhyperbolic fixed points that might be 
late-time attractors (for instance points having at least one unstable eigen-direction, like Q22) 
are excluded from the analysis). These are the expanding solutions Q 20 and Q 21 and the 
contracting one R^ 2 - We remind that the points Qi and their contracting partners points 
Ri through the symmetry (3.50), exhibit opposite dynamical behaviors, and hence from the 
following analysis we also obtain information for the contracting solutions -R20,-R21 and the 
expanding one <522- 

In order to calculate the center manifold of Q 20 = (1,0,1) for the system (3.48) we 
introduce the variables 

e=l — Ti, u = l-Q 0 , v = Q de , (A.57) 
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which satisfy the evolution equations 


(A.58a) 
(A.58b) 
(A.58c) 


-3(1 - u)( 1 - e)e 2 , 

-^(2 - u)u [e(3 v(w m + 3) - 3w m - 1) - 6u] , 

3(1 -u)( 1 - v)v [(w m + 3)e - 2] . 

The center subspace of the origin of (A.58) is spanned by the vectors (1,0,0) and (0,1,0), 
which implies that the local center manifold of the origin can be written locally as the graph 
{(e,u,v) : v = h(e,u),h( 0,0) = 0,Dh(0,0) = 0 , ||(e,it)|| < 5}, with 6 a suitably small 
constant and Dh the matrix of derivatives. The function h(e, u) satisfies the quasilinear 
partial differential equation 

- ]-(u - 2 )u^- [3 {(w m + 3)e - 2} h - ( 3w m + l)e] + 3(tt - l)(e - l)e 2 ^ 

2 ou oe 

+ 3(u - 1) {( w m + 3)e - 2} (h - 1 )h = 0, (A.59) 

which admits the formal general solution 


(A.60) 


(-In 

(u — 2 )u(l — e ) 2//3 

,ln 

e _2 / e (e - 1)(1 - e)^-™-- 1 (h(e,u) - 1 ) 

V 2 " 

e 2 / 3 /i(e, u) 

1 

rb 

1_ 


de 

df 

du 

df 

dv 

df 


and the trivial solution h = 0. Nevertheless, in order to complete the analysis numerical 
investigation is required. Using Taylor expansion we obtain that, up to third order in the 
vector norm, the solution of (A.59) is the trivial solution h = 0. Thus, the dynamics on the 
center manifold is given up to the same order by 

=-3{1-u)e 2 + 0(3), (A.61a) 

df 

du 

— = u(3w m + l)e + 0(3). (A.61b) 

dr 

Neglecting the error terms, and eliminating the time variable, we obtain the equation 


= 3(u - l)e(u) 
1 j u(3w m + l) 


(A.62) 


which admits the solution 

3(u — ln('Li)) 

e(u) = c\e 3™m+i , (A.63) 

that satisfy lmp^o = cioo, and hence it is infinity unless c\ = 0. Therefore, Q 20 is a saddle 
for w m > 0, as shown in Figure 6. Finally, by symmetry, f?20 behaves as saddle too. 

In order to examine the stability of the curve of critical points Q 21, using the center 
manifold theorem, we introduce the variables 


u\ — Qoc Q 0 T 


d 1 - Qoc ) ®DE 

2Qoc 


,u 2 = - (1 - Qoc ) (1 - 7i)(l + 3 w m ),v = Ode, (A.64) 
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where Q q c €E (0,1) is a constant, which satisfy the evolution equations 


du\ 

df 

+ 

dU 2 

df 

dv 

df 


3v(v - 1) [Ql c (v - 2) + 2Q 0c m - v ] [w m (3 Ql c + u 2 - 3) + Ql c + 3u 2 - l] 


2Q()c{3'w rn + 1 ) 


QocV . V . 

— + 2 Q^ + Qoc - Ul] ~ 1 


u 2 [-3v(w m + 3) + 3w m + 1] 

(Qlc ~ -*-) (3 w m + 1) 

3u\ [Qoc( v - 2 ) + 2Qoc«l - v] [3 (Qoc - 1) Wm + Q\ c + 2 U 2 - l] 

Qoc ( Qlc - l) 2 (3rc m + l) 2 

O , , , ( QocV , v 

-3v(v - 1 )---h 777 -h Qoc - Ui 


— 3v 


2 Qoc 


2 u 2 (w m + 3) 


(Qlc l) Q> w m + 1) 


- 2 


(A.65a) 
(A.65b) 
(A.65c) 


The center subspace of the origin of (A.65) is spanned by the vectors (1,0,0) and 
(0,1,0), and hence the local center manifold of the origin can be written locally as the graph 
{(ui,u 2 ,v) : v = h(u\,u 2 ), M0, 0) = 0, Dh(0, 0) = 0 , ||(ni, u 2 )\\ < <5}, with 5 a suitably small 
constant, and where the function h satisfies the quasilinear partial differential equation 


3ul [3 (Ql c -l)w m + Qlc + 2 u 2 - l] [(Ql c - 1 )h + 2 Qq c {u\ - Q 0c )] dh 
Qoc(Q 2 0 c -l) 2 (3w m + l) 2 du 2 

f 3 [w m (3 Ql c + u 2 - 3) + Ql c + 3 u 2 - l] (h - 1 )h [{Ql c -l)h + 2Q 0c (m - Q 0c )] 


+ 


2 Qoc 

2 u 2 (w m + 3) 


Qoc + 


- 1 


(Qoc l) Q> w m + 1) 


+ 2 


(h — 1 )h 


2Ql c {3w m + 1) 

u 2 [-3 (w m + 3 )h + 3 w m + 1] 
{Qo c ~ l) Q> w m + 1) 

(Qoc-i)^ 


— 3 h 


dh 

du\ 


2 Qqc 


T Qoc 


= 0. 


(A.66) 


Assuming that h(u\,u 2 ) = anu 2 + 012^1^2 + 022^2 + 0(3) and evaluating up to third order, 
we find a good approximation of the center manifold given by 

h(ui,u 2 ) = -3a 12 Q 0c u\ + ai 2 uiu 2 - 777^2 + 0(3). (A.67) 

uCtf 0c 


Therefore, neglecting the 0(3)-terms, the evolution on the center manifold is governed by 
the equations 


du\ 

2Qoc^1^2 

(A.68a) 

df 

U2 1 ~ Qlc ’ 

du 2 

QQocul 

(A.68b) 

df 

(l — Qq c ) (3w m + 1) 

where by definition, u 2 > 0. 




Integrating (A.68) we acquire 


2c 2 Qo e [ci ( Ql c - l) (3 w m + 1) + 6Q 0c t] Wrn+ l + Q\ c - 1 

2Qoc 

_ (Qoc ~ -*-) Q> w m + 1) _ 

(Qlc - !) ( 3w m + 1) + 6 Q 0c f 


(A.69a) 
(A.69b) 
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Taking the limit f —»■ oo in the above expressions we obtain (u±, u 2 ) —>• (C200,0), C2 / 0, u 2 7^ 
0. In the special case where C2 = 0 we obtain the limits ( 111 , 112 ) —> ( °q q , 0)■ In both cases 
the origin is unstable along the iq-axis. Since it is stable along the U2-axis, it follows that 
Q 21 is a saddle. 

Another nonhyperbolic point that can be a late-time attractor is the contracting point 
R 22 , which has a 2D stable manifold. Introducing the new variables 

e = 1 - T\,x = Qo + l,y = 1 — Qde, (A.70) 

we obtain the equivalent dynamical system 

— = 3(x — l)(e — l)e 2 , (A.71a) 

df 

dx 1 

— = -(x- 2)x{3y[(w m + 3)e - 2] - 8e + 6}, (A.71b) 

ctr Z 

= 3(x - 1 )(y - 1 )y[(w m + 3)e - 2], (A.71c) 

where the local center manifold of the origin (e, x, y) = (0, 0,0) is tangent to the e-axis. 
Hence, it can be written locally as the graph 

{(e,x,y) : x = hi(e),y = h 2 (e),hi(0) = 0,h 2 (0) = 0, h! x (0) = 0,/i' 2 (0) = 0, |e| < <5}, (A.72) 


where <5 is a suitably small number. The functions h\ and h 2 must satisfy the quasilinear 
system of differential equations 

(ft! - 2)h\ {3 h 2 [(w m + 3)e - 2] - 8e + 6} - 6(e - l)e 2 (hi - 1 )h[ = 0, (A.73a) 

(hi - 1) { (h 2 - l)h 2 [(w m + 3)e - 2] - (e - 1 )e 2 h x } = 0. (A.73b) 

These equations admit the following classes of solutions: 

1. The trivial solution hi = 0, h 2 = 0, 


2. the one-parameter class of solutions 


hi 



h 2 — 0, 


e 2/3 e 2( Cl + i) 

(1 — e) 2 / 3 ’ 


(A. 74a) 
(A.74b) 


3. the one-parameter class of solutions 

hi = 0, (A.75a) 

h '2 = -2---, (A.75b) 

e Cl ~~(l — e) Wrn+1 e~ Wrn ~ 1 + 1 


4. the 2-parameter class of solutions 


hi 

h 2 


1 ± 


I g2c 2e 2/3(l _ g^m+g 

e Cl (1 — e) w m + 1 _|_ g2/ e^w m +1 

1 


+ 1 , 


g Cl_ 7 (1 _ gitUm+lg—TO m —1 _|_ l 


(A.76a) 
(A.76b) 
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From all the above solutions the only one that satisfies the conditions /ii(0) = 0, /i2(0) = 
0, h\ (0) = 0, h' 2 (0) = 0 is the trivial one. Thus, the dynamics on the center manifold is given 
by 

^ = 3(l-e)e 2 , (A.77) 


that corresponds to a gradient-like equation for the potential 17(e) = l/4e 3 (—4+3e) for which 
the origin is a local maximum. Thus, the center manifold of the origin is unstable, and hence 
R 22 is a saddle. Integrating out the above equation we extract the solution 


f (e) = ci - \ tanh '(1 - 2e), (A.78) 

3e 3 

for which the origin is approached as f —>• —00 and not as f —» +00, that indeed confirms 
the above statements. 
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